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Abstract. In order to describe forward hadron productions in high-energy nuclear 

. collisions, we propose a Monte-Carlo implementation of Dumitru-Hayashigaki- Jalilian- 

f->i ; Marian formula with the unintegrated gluon distribution obtained numerically from 

' I the running-coupling BK equation. We discuss influence of initial conditions for the 

flj ■ BK equation by comparing a model constrained by global fit of small-x HERA data 

,S^ , and a newly proposed one from the "running coupling MV model." 

T-H . 

^ ! 1. Phenomenology with rcBK equation 

m 

en ■ Relativistic Heavy-Ion Collider (RHIC) at Brookhaven, the QGP machine, has provided 

a new opportunity to explore the physics of the dense gluonic state at small Bjorken's 

Q ' X in incoming hadrons, which is called Color Glass Condensate (CGC) [1]. Analyses 

of (e.g.) deeply inelastic scatterings (DIS) at HERA[2, 3] and heavy-ion collisions at 
RHIC [4] support the phenomenology based on the CGC picture. 

L^ ■ A recent highlight is the inclusion of the running coupling effects in the Balitsky- 

^ • Kovchegov equation (the rcBK equation) [5, 6], which controls the x-evolution of the 

unintegrated gluon distribution J\f{k, y) in the incoming hadron {k is the transverse 
momentum and y = ln(l/x) is the rapidity), which is obtained as the Fourier-transform 
of the 'dipole scattering amplitude' J\f{r, y) with r being the transverse size of a dipole. 
The rcBK equation has been confronted with the global analysis of DIS [7] as well as with 
the data obtained in proton-proton (pp), proton- nucleus (pA) and nucleus-nucleus (AA) 
collisions [8, 9, 10, 11] (see [12] for recent progress). Our goal is to elaborate a detailed 
model for particle productions in forward regions at RHIC and the Large Hadron Collider 
(LHC), by applying the gluon distribution obtained from the rcBK equation. 

2. Forward particle production 

2. 1 . Initial conditions for rcBK equation 

With the DIS data for x < xq accumulated at HERA, one can constrain the numerical 
solution of the rcBK equation quite accurately. This program was accomplished in [7] 
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by choosing an initial condition at x = Xq = 0.01 as (e.g.) 



A/'(r, yo) = 1- exp 



4 I Ar 



(1) 



Among several parameter sets now available, we adopt here the AAMQS parameter set 
h {Q% = 0.1597 GeV2, 7 = 1.118, C = 2.47, A = 0.241 GeV, afr=l) [7]. 

In addition, we will try a new initial condition, motivated by the McLerran- 
Venugopalan model with the running coupling modification (rcMV) [13]: 



A/'(r, yo) 



afi 



1 — exp 
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where a,(r) = l/[boln{4:C^/{r^A^) + e)] with 60 = {ll-2N{/3)/4:n with iVf = 3, and a^ 
is fixed by requiring that maximum value of A/'(r, i/q) be unity. We have set Af{r, hq) = 1 
for r larger than this maximum point. For the rcMV initial condition we use the 
parameters C = 1, A = 0.2 GeV, Q\ = Q^J \n{QlJA'^) with Q^^ = 0.2 GeV^ as a trial. 

2.2. DH J formula 

In particle productions at very forward rapidities y > one can probe the saturation 
regime of one of the colliding hadrons. In this situation a factorized formula for the 
cross section is proposed by Dumitru, Hayashigaki, and Jalilian-Marian (DHJ)[14]: 

dN K ^^ [^ dz 



dyhd'^PT (27r 



|2 A^ 



pi 1 

/ — xifi/p{xi,p'^) Sf ( — ,X2) Dh/i{z,pl) , (3) 

Jxp ^ \ Z / 



where fi/p is the collinear distribution function for the large-Xi parton i, M describes the 
small- 0:2 unintegrated gluon distribution, and Dhji deals with the high-p^ fragmentation 
of a parton i into a hadron h with the momentum fraction z. The parameter K is 
introduced to absorb the higher-order contributions. In AA collisions, /j/p and M are 
to be generalized to those for nuclei. 

The formula (3) was applied [14] to forward hadron productions in d-Au collisions 
at RHIC, with adopting a certain CGC model for M . The larger gluon density in 
a nucleus is advantageous for probing the saturation effects. Recently Albacete and 
Marquet[9] exploited the rcBK equation to improve the M part theoretically. Assuming 
a homogeneous nucleus in the transverse plane, they obtained the best description to 
date for the rapidity and momentum dependence of the d-Au data at RHIC. 

2.3. Monte-Carlo implementation for nuclear collisions: the MC-DHJ/rcBK model 

For more quantitative study, we need a nuclear model including the impact-parameter 
dependence and initial fluctuations. In a MC implementation developed in [15], nucleons 
are randomly sampled according to the Woods-Saxon density in an event-by-event basis, 
and the initial saturation momentum Q^q ii^ (1) ^'^'^ (2) is taken to be Q^q = N{ri_)Qlp at 
each transverse coordinate, where N{r±) is the number of nucleons in the dense target 
within a transverse area S^. We apply locally in the transverse plane the numerical 
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Figure 1. Comparison of forward hadron spectra computed with the AAMQS set h 
(solid) and rcMV (dashed) to the data observed at RHIC. The K factor is chosen as 
K — 1.5 for charged particles (/i~) and K = 0.5 for tt^ in pp and d-Au collisions. 



solutions for the rcBK equation with different initial values of Q^q. Such an approach 
has been pursued within the so-called k^ factorization approximation for AA collisions, 
which reproduces successfully the centrality dependence of the hadron multiplicity at 
RHIC and the LHC[15, 10, 11]. 

Here we combine the DHJ formula with the rcBK evolution in the MC 
implementation (MC-DHJ/rcBK). That is, we compute particle productions at each 
transverse grid r^ using the DHJ formula (3) with A/'(A;, y) numerically obtained from 
QIq at grid r^ determined by the MC code: 



dN 



Td{ri_) 



dN 



X 



(4) 



DHJr± 



dyhd'^PTdr± dy^d'^px 

Here Td{r±) is the thickness function on the dilute side. We stress as an advantage of 
this approach that there is no more additional parameter after fitting pp collisions. We 
comment also that the MC implementation allows us to study the initial fiuctuations[15]. 



3. Results 

We use for fi/p and D^/i the CTEQ6M NLO PDF [16] and DSS NLO fragmentation 
functions [17], respectively, and set the factorization scale to /x^ = p^. We remark that 
an oscillation appears in N{k, y) for smaller QIq when a sharp cutoff for the running 
coupling asir) = l/[6oln(4C^/r^A^)] is adopted at a^. Thus we tried a smooth cutoff 
as{r) = l/[6oln(4C^/r^A^) + a] where constant a is adjusted to make as{r) — )• 2 as 
r — )■ oo in the rcBK evolution in the case of the rcMV initial condition. 

In figure 1, transverse momentum distributions of negatively charged hadrons 
h~ at pseudo-rapidities 7] = 2.2 and 3.2 from BRAHMS [18] and neutral pions 
TT^ at rj = A from STAR [19] in pp and d-Au collisions at ^/s = 200 GeV 
are compared to our results. The AAMQS set h with K = 1.5 (0.5) describes 
the forward particle multiplicities of h~ (vr^) very nicely in pp and d-Au collisions 



Forward particle productions at RHIC and the LHC from CGC with local rcBK 



at the same time without changing any parameters 
hand, the rcMV initial condition with our current 
(0.5) for h" (7r°) leads to a good agreement with 
but not for pp collisions. It is very 
interesting to do global analysis of HERA 
data by using rcMV initial condition and 
fix the parameters. We will report more > 
systematic analyses on the parameter j: 
dependence elsewhere. f 

Finally a result from test calculations 
of the MC-DHJ/rcBK at extremely for- 
ward rapidity r] = 8.5 for pp collisions at 
a/s = 7 TeV, which is measured in LHCf 
experiment [20], is shown in figure 2. One 
expects that dependence of the rcBK evo- 
lution on initial conditions will become 
weaker at higher energies, but we see that 
the result is still relatively sensitive to the 
initial conditions. 



in the model. On the other 

parameter set with K = 1.5 

the data for d-Au collisions, 
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Figure 2. Transverse momentum 
spectra of average of charged pions 
(bold) and of charged kaons (thin) 
at 77 — 8.5 for pp colhsions at 
yi = 7 TeV from MC-DHJ/rcBK 
model (^=1) with the AAMQS set 
h (sohd) and rcMV (dashed). 
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